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Abstract 

We reconsider the Euler-Lagrange equation for the Skyrme model in the hedgehog ansatz and 
study the analytical properties of the solitonic solution. In view of the lack of a closed form solution 
to the problem, we work on approximate analytical solutions. 

We show that Pade approximants are well suited to continue analytically the asymptotic rep- 
resentation obtained in terms of a power series expansion near the origin, obtaining explicit ap- 
proximate solutions for the Skyrme equations. We improve the approximations by applying the 
2-point Pade approximant procedure whereby the exact behaviour at spatial infinity is incorpo- 
rated. An even better convergence to the exact solution is obtained by introducing a modified 
form for the approximants. The new representations share the same analytical properties with the 
exact solution at both small and large values of the radial variable r. 
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I. INTRODUCTION 



The Skyrme model provides a picture of baryons interacting via meson exchanges p| , 
0. The model is based on the pre-QCD nonlinear a model, with the same group of chiral 
symmetry G = SU{Nf) x SU{Nf) [Q. Classical stability arguments require the presence 
of an additional term in the nonlinear cr-model Lagrangian this term was introduced 
initially by Skyrme. In this scenario, baryons emerge as solitonic solutions of the Euler- 
Lagrange equations which are highly nonlinear second order equations and cannot be 
solved in exact closed form. Solutions are obtained by numerical integration. This is an 
awkward situation for the evaluation of physical quantitie s. In practice it is more convenient 
to have explicit approximate representations for the solitonic solutions. 

In the present work we obtain reliable approximate representations for the solution of the 
Skyrme model in the SU{2)l x SU{2)r sector by means of the Pade Approximant (PA) 
procedure whereby the series expansion solution of the differential equations is continued 
analytically. 

The successful implementation of PA summation relies on the simple structure of the 
series solution near the origin which is almost an alternate geometric series. In order to 
reproduce the essential features of the exact solution at large distances, we use the 2-point 
PA method and construct modified representations that share the exact properties of the 
skyrmion configuration near both boundaries. Futhermore our approximate solutions are 
reliable in the whole interval of the independent variable r, namely r = to infinity. 

This paper is organized as follows. In section II we introduce the Skyrme model in the 
SU{2)l X SU{2)ii sector and restrict the problem using the hedgehog ansatz. In section III 
we present properties of the hedgehog solution showing power series r epresentations for the 
Chiral Angle. In section IV we introduce and use our approach based on the PA method. 
Finally, in section V we draw the conclusions of our work. 

II. THE SKYRME MODEL 

Following Skyrme's proposal, baryons are soliton in the cr— nonlinear model with an 
additional stabilizing term. We will restrict ourselves to the sector SU{2)l x SU{2)ji where 
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the Lagrangian density of the model reads 
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(1) 



being U a unitary operator. The parameter /,r is the usual pion decay constant whose 
experimental value is 93 MeV and g is the dimensionless Skyrme parameter. Equivalently, 
U can be expressed in terms of a scalar and a pseudoscalar field, a and tt respectively, in 
the following form 

U={a + tr-7z)/U (2) 

It should be stressed that a and tt are not independent fields, and are related through 
0"^ + TT^ = This restriction is the source of nonlinearities in the theory. 
Rewriting tt and a in terms of a new field F = fnFip 



cr = cos F, 77 = f^if sin F, 



and going back to U, we obtain the representation 



U = exp{iT ■ F) 



(3) 



(4) 



Previous analysis ^ have shown that a reasonable agreement with experimental values of 
physical quantities is obtained when the hedgehog configuration of U is adopted, this is a 
particlular static ansatz where F takes a spherically symmetric form, namely 



Us{r) = exp(iT • n F{r)) 



(5) 



where h denotes a unit radial vector. 

In this representation, the energy reduces to a functional of the Chiral Angle, F{r), alone 



E[F] = 271 df fy (^^^ ' + 2fl sin^ F + 1 sin^ F L 



'dF\ sin^F 
df I 



(6) 



The variation 6E/SF = generates a nonlinear differential equation for F{r), whose 
solution corresponds to the static configuration of minimal energy. Written in terms of the 
dimensionless variable r = gf-j^r, the differential equation for F{r) reads 



[r^ + 8 sin^ F{r] 



d^Fjr) 
' dr^ 



— sin^ F(r) sin[2F(r)] + sin[2F(r)] - 2r 



dF{r) 
dr 



-4 sin[2F(r)] 



dF{r) 
dr 



1 2 



(7) 
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As finite energy solutions are required, U{r) must tend to an arbitrary constant element of 
SU (2) at spatial infinity. Choosing f/(r) ^ 1 as r ^ oo , implies the boundary condition 

F(oo) = (8) 

This condition defines a mapping S^, from the compactified configuration space 

to the identity in the target space in SU{2) which is isomorfic to 5''^, falling into distinct 
equivalence classes labeled each by the winding number Z 

MSUi2)) = U;{S^) = Z, (9) 

The integer number Z, that counts the coverings of the target space, is defined as the 
topological charge q = J (Ix^Bq related to the baryonic current 

B, = ^^e,^^^Tr[{U^d''U){U^d^U){U^d'^U)] (10) 

For the hedgehog ansatz, the baryon number charge density is 

1 dFsin'^F 

^°-~2^^^^ ^^^^ 
and corresponds to a baryon number 

B = — [2F(0) - 2F(oo) - sin 2F(0) + sin 2F(oo)] (12) 

From the foregoing, it follows that the solution belonging to the sector with baryon number 
B = 1, satisfies the condition 

F(0) = TT. (13) 

The numerical solution labeled with i? = 1 is shown in Fig.l . Inserting this solution into 
Eq. (D yields E = 23.27TU/g. 



III. PROPERTIES OF THE HEDGEHOG SOLUTION 

Before applying the PA procedure to the Skyrme model, we shall consider the solution 
of the Euler-Lagrange equation at both ends of the definition domain. Near the origin, a 
power series solution can be obtained 

Snir) =n + F,r + ^F^r^ + ^F.f' + ... , (14) 
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where Fi turns out to be undetermined. All even powers have vanishing coefficients. Odd 
power coefficients are written in terms of Fi and the dimensionless parameter 0, defined as 
= Fi/gfj,] for example, 

^ 24 1 + (32/5)0^ + (88/5)0^ + (488/5)0^ 
^ 7 ^ 1 + 2402 + 19204 + 51206 ' ^ 

while the results for higher terms may be found elsewhere |p. In terms of the dimensionless 

variable r = gfT^f, expansion (14) becomes 



^^(r) = 7r + 0r + ^03r3 + i05r^ + ..., (17) 

where 03 and 05 are given by the same expressions as F3 and F5 but with instead of Fi. 
From now on, we shall be concerned only with the dimensionless variables. 

For large values of r , there is a solution of the form C/r^, as it can be readily seen from 
trying this form in equation (|^) and keeping only linear terms in F{r). 

In order to find contributions of higher order in 1/r, it is appropriate to perform the 
change of variable p = 1/r in equation (|^. In this new variable it reads: 

[p2 + 8pSinF(p)]^^^|^ + 8p='sin2F(p)^ + 4pSin2[2F(p)] ^"^^^^^ 



dp^ dp 



dp 



-Ap^ sin^ F{p) sin[2F(p)] - sin[2Fp] = (18) 

Trying a power series expansion as before, 

F(p) =Ko + K,p + ^K2p' + ^K^p' + ... (19) 

we find that the series solution contains just even powers. The first few coefficients written 
in terms of K2 are 

= (20) 

K, = -^-^Kl (21) 

Ks = -6720KI (22) 

Just as the initial slope 0, the parameter K2 cannot be determined from the differential 
equation. Approximate values for and K2 are given in Ref. [^]. These values were deter- 
mined through the use of integration routines based on shooting methods H] . This is some 



sort of eigenvalue problem for and K2. The solution is obtained by numerical integration 
of the differential equation, starting from the origin where the value of (j) must be such as to 
satisfy the correct boundary condition corresponding to the relevant solitonic solution. Once 
(f) is found, the value of K2 is obtained in similar way using equation (|TH). The reported 
values are . 



= -1.0037 K2 = 17.2772 (23) 

These values, as well as the vality of the numerical solution can be checked in a consistent 



manner. We will not go on through this and just refer the reader to the Ref. |jTT 

Let us make a slight modification in our analysis by rewritting the power expansion (|1 
in the following form 



S^ir) = TT + 0r + air^ ll + ^r' + ^^r' + ...), (24) 

^ ai ai 02 ' 



where the coefficients ai, a2, 03, ... are defined through 

ai = ^, "2 = ^ as = ^, (25) 

and so on. As shown in table I, the numerical values of the ratios — a„+i/a„ vary very slowly 
with n. Thus, one can introduce the nearly constant R ~ —an+i/cLn (for any value of n). 



Consequently, the portion within the brackets in (^) can be replaced, to a good degree of 
approximation, by the alternate geometric series 

1 - i?r2 + {Rr^ - {Rr^f + {Rr^f - (26) 

which has a radius of convergence tq defined by the relation Tq = |1/-R|. 



Within the radius of convergence, r < ro, expansion (|2^) may be replaced by 



~ TT + 0r + air\l - Rr^ + {Rr^ - {Rr^ + [Rr^y - ...). (27) 

An alternate series like l — x + x'^ — x^ + ... can be continued beyond its radius of convergence 
by the function (1 + x)^^, so that Sh can also be analytically continued for all positives 
values in the real r axis using the representation 

F{r) = vr + 0r + (28) 

Taking R = —02/01, Eq. (pSD is exactly the [3, 2] order Fade Approximant representation of 
the series {^^- The obvious task is to test whether higher order of approximation can provide 



a better representation of the soliton. There is a great amount of examples [|I2| where Pade 
Approximant procedure provides rehable analytic continuation of power series expansions 
and that PA representations converge rapidly to the relevant function of the problem. In the 
next section we will build different representations of PA to the hedgehog solution showing 
that the features of reliability and fast convergence also hold in this problem. 

IV. FADE APPROXIMANTS TO THE HEDGEHOG SOLUTION 

Our approach is based on Pade approximants, which are rational functions used to provide 
an analytic continuation of a power series representation of a given function, and are typically 
known to accelerate the convergence of the series. The Pade approximant P[m,n]{x) of order 
[M, N] to the series S{x) = J2n dn^" is defined as the ratio of two polynomials 

M 

P[mM^) = ¥ (29) 

where we set Bq = I without loss of generality. The remaining M + N + 1 coefficients are 
chosen so that the first M + N+1 coefficients in the Taylor expansion of P[m,n] coincide with 
the series S{x) through order M + N. Conversely, if only the first M + + 1 coefficients 
of the series S{x) are known, such a Pade approximant can be used to predict the next 
coefficient in the series IjlSj . 

The PA of order [M, A^] to the hedgehog solution is found by comparing its Taylor ex- 
pansion at the origin with series ([T^) . From all the possible combinations {M,N) allowed 



by the input series expansion (JT^, we are interested in tho se that may provide a suitable 
representation to the hedgehog solution. In order to enforce the conditions that fix the 
soliton solution to be of winding number one, we will require that the behaviour at infinity 
of the PA representations is given by the leading term of expansion (p^Qf) . To this end we 
impose the contraint N — M = 2. We have computed the values of coefficients a/c and bk 
corresponding to the sequence fo PA's [0,2], [1,3] and [2,4] using the value 0i = —1.003 
found by numerical methods. The results are displayed in figure 1. The explicit dependence 
of the PAs in terms of the coefficients of series (|T^ is 

TT 

^[0,21 = : 7—72' 



for [0, 2] and [1, 3] respectively. 

The PA representations built from the exact series solution near the origin do not guar- 
antee the exact behaviour at intermediate and large values of r. However, one could use the 
asymptotic solution at infinity given in (^) in order to find better representations of the 
hedgehog [Q. Yet, the right behaviour at infinity may be incorporated using the 2-point 
Pade Approximant procedure. This is a natural extension of the usual PA method which 
consists in fixing almost one coefficient of the approximant using the asymptotic series solu- 
tion at infinity 0. We have computed the sequence of 2-point PA of order [1, 3], [2, 4], [3, 5] 
(see figure 2). Certainly, the 2-point PA representations are more reliable than the ordinary 
ones. The explicit form of PA [1, 3] is 



7r 



J<:2(6gi3+</,3^2) 
6(ii'2</>2-27r3) 



TP _ <oyi\2<p--^T^-) /o-|\ 

_ 7r(iC2(/>3 + 12gi7r) _ (/'(i<:203 + 1207r) , 6.^3+037r^ ^ ' 

The next step is to include the subleading behaviour of the exact solution at infinity 
within the PA representations. We have done this up to order 0{l/r^) as shown in figure 
3. As we go on with the procedure, the degree of approximation is improved both near the 
origin and for large values of r. However, for intermediate values of r, the PA representations 
are slowly convergent. Nevertheless we will show that there is a suitable functional form 
that provides very good representations to the hedgehog solution. 

To this end, let us prescribe a modified PA of order [0, 4] through 

^[0,4] + + ^2^2 + 5g^3 + ' (^2) 

such that the coefficients oq, &i, &2 and 63 are determined by expanding F[o,4](r) in a power 
series near the origin and matching the terms, order by order, with those of the exact 
solution (p!7|) up to O(r^), and the remaining coefficient 64 is fixed by using the leading term 
at infinity given by (p^Qf). We found explicitly, 

71 

" (I 2</,^ . 3<p2^2 I -12</.3-fe7r2 ^3 . 4,r2 ^4^1/2" ^^^^ 

TT TT^ StT-^ K2 ' 

The form (|3^ ) can be generalized to higher orders [0, k] as 

^f^''^^""^ " {l + hir + h2r^ + ... + hkr^f/^ ^^^^ 



where the unknown coefficients may be found following the prescription described for the 
particular case k = 4. We have built the representations corresponding to = 6, = 8, 
and k = 10 (see Fig. 4). By construction, such representations have the exact Taylor 
expansion near the origin up to order O(r^), O(r^) and O(r^) respectively, as well as the 
exact leading behaviour at infinity. The new PA representations improve in a remarkable 
way the representations obtained so far. 

Following the program, an even better agreement for the large r behaviour may be ob- 
tained by including in the modified representations (|3^) the subleading behaviour of the 
hedgehog solution. In figure 5 we have depicted the approximants so constr ucted. The 
figure shows an exceptional good agreement betwen the numerical solution and the modified 
PA given by (0). In order to show the utility of PA procedure as a method of analytic 
continuation, we have included in the same figure the power series representation near the 
origin up to and including order O(r^). 

The utility of PA method in this problem is outstanding, particularly for its simple and 
reliable application to the highly nonlinear Euler-Lagrange equation of Skyrme. Moreover, 
the PA representations reproduce exactly the main properties of the hedgehog solution at 
both small and large values of r. 

In order to further check the reliability of the PA representations, we have calculated its 
corresponding baryon number given by 

fi = -— / —Anr'dr, 35 

ZTT^ JO dr 

and the minimal energy (P) whose numerical value is E = 72.88f-^/g. We obtain B = 1 
for the 2-point PAs which is exactly the topological number of the skyrmion. As for the 
energy, the results for the last sequence of approximants in units of /t^/q are -E[o,4] = 73.60, 
Eiofi] = 73.02, -^[0,8] = 72.96, and -E[o,io] = 72.94, which converge to the numerical value 
reported in previous work s [Q. 



V. FINAL REMARKS 



We have shown that the suggested PA representations for the Skyrme solutions can 
incorporate in a simple way the properties of the exact solution near both physical boundaries 



of the problem. This solves the difficulty acquainted in Ref. [11|, wh ere the authors assert 
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the impossibility to obtain a PA to the series solution about the origin which, in addition, 
reproduces the exact behaviour of the chiral angle at infinity. Moreover, our representations 
avoid the use of two separate analytical repr esentations for the skyrmion, one for each of 
the two regions, which was the common approach to this problem |T^, clearly unsuitable 
for phenomenological calculations. 

Certainly, the results presented here highlight the utility of rational functions, PA in 
particular, for analytic continuation of asymptotic series. In this problem, the knowledge of 
the asymptotic behaviour at two points favors a better convergence to the relevant function. 
The approximate analytical solutions can be used with reliability in the evaluation of phys- 
ical quantities. The calculations are further easier and more manageable than numerical 
computations. 
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TABLE I: Ratios of the coefficients of the series expansion 



02/01 03/02 04/03 05/04 oe/os 07/06 03/07 
-0.102 -0.136 -0.133 -0.139 -0.147 -0.149 -0.152 
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FIG. 3: Two point PA exact up to order 0(l/r^) in the limit r oo. 
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FIG. 4: Sequence of modified PA to the hedgehog solution. 
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FIG. 5: Modified PA representations which reproduce the exact leading and subleading behaviour 
at r — > oo. 
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